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We describe application of the gauge/gravity to study of depairing current in thin superconducting 
wires. The large number N of colors of the gauge theory is identified with the number of filled 
transverse channels in the wire. On the gravity side, the physics is described by a system of D3 
and D5 branes intersecting over a line. We consider the ground state of the system at fixed electric 
current and find that there is a continuous phase transition at a critical current, in the universality 
class of the dissipative XY model. We discuss relation of our results to recent experiments on 
statistics of the switching current in nanowires. 



Destruction of superconductivity by a current flowing 
through the superconductor has long been recognized as 
an important topic both from the viewpoint of funda- 
mental physics and for applications (for a review, see 
Bardeen One expects that theoretical treatment 

should be the simplest for samples that are effectively 
one-dimensional (Id), i.e., wires in which the supercon- 
ducting density depends on only one coordinate. Even 
for this case, however, a complete theoretical treatment 
of the transition at fixed current has not been forthcom- 
ing. The main difficulty lies, somewhat ironically, in the 
definition of a current-carrying normal state: the equilib- 
rium Fermi distribution used by the conventional mean- 
field theory is clearly inadequate for the purpose as it 
carries no current at all. 

A related but different problem is a transition caused 
not by a fixed current (sustained by an external battery) 
but by an initial winding of the order parameter. Ex- 
perimentally, this is the condition appropriate for a thin 
superconducting ring. In this case, the supercurrent is 
only metastable: there are fluctuations that lower the 
free energy [2| (provided there is an amount of disorder 
or a finite temperature Q ) ; these have become known as 
phase slips. The relevant question to ask then is whether 
there is a maximal winding number density beyond which 
the metastability becomes classical instability. 

To define a current-carrying normal state, one needs to 
include some mechanism that equilibrates the electrons 
with respect to momenta. One possible approach is to 
include weak scattering via a kinetic equation. Another, 
which we adopt here, is to start at the opposite extreme — 
a theory with strong electron-electron interactions. One 
may hope that for such a theory there is a complementary 
(dual) description in terms of weakly coupled collective 
modes. If these collective modes are the same as seen 
on the superconducting side, one will have a unified de- 
scription applicable to both phases, which should make 
understanding the phase transition easier. 

Recently, following the discovery of the gauge/gravity 
duality carrying out this program has become prac- 
tical. The gauge/gravity duality allows one to study a 
strongly coupled SU (N) gauge theory with a large num- 
ber N of colors by doing calculations in classical gravity, 



albeit in a higher-dimensional spacetime. Here we use 
this method to study the depairing transition at fixed 
current. The number N of colors is taken to correspond 
to the number of populated transverse modes (channels) 
in the wire. Our conclusion is that at a critical current 
there is a continuous phase transition to the normal state 
in the universality class of the dissipative XY model. The 
duality can also be used to study transition due to an ini- 
tial winding; we outline the strategy in what follows. 

In the construction we use here, the SU(N) gauge 
fields live in (3+1) dimensions, while the electrons only 
move in (1+1). The (3+l)-dimensional sector is confor- 
mal: qualitatively, this corresponds to an attractive 1/r 
potential between left- and right-moving electrons. Re- 
gardless of whether one considers this a realistic model 
of a Id superconductor, the main result — a second-order 
transition at fixed current — may be universal enough to 
apply to experimental samples. 

The large N suppresses interaction of the collective 
modes at the cutoff scale; in a wire, this corresponds 
to suppression of self-interaction of the pair field by the 
wire's cross-sectional area. As a result, in a given experi- 
ment, the near-critical behavior may be controlled by the 
ultraviolet Gaussian fixed point, rather than the stable 
infrared fixed point of the dissipative XY model. As we 
discuss at the end, this appears to be the case in recent 
experiments 0-0]' a ^ least at not too low temperatures. 

We begin by establishing our convention for assembling 
electron operators into Dirac spinors. Assuming that su- 
perconductivity is due to a correlation between oppo- 
sitely moving electrons in the same transverse mode, we 
expect it to show in correlation functions of the opera- 
tor Ea=i a R^L^ where the subscripts R and L designate 
the electron operators with positive and negative mo- 
menta, respectively. For this to be a color singlet, 
should transform as N of the SU(N) and &z, as N (or 
vice versa). Thus, we define a 2-component Dirac spinor 
yj as follows (omitting the color index A): 



(1) 



This is in the representation where the Dirac 7 matrices 
are given by 7 = ax, 7 1 = — 102, and 7 s =173, in terms 
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of the Pauli matrices <x. The identification of the super- 
conducting channel as apbL implies that a} Rk creates a 
k > electron, and b\ k a k < electron (i.e., o) Lk creates 
a k < hole). With this convention, the upper and lower 
components of tp have opposite electric charges, the su- 
perconducting channel is (where tjj — ip^j ), and an 
external electromagnetic potential couples to the axial 
current ^7 M 7 5 ?/>. 

To reproduce the physics of such a superconductor on 
the gravity side, we consider the system of N coincident 
D3 branes and a single D5 brane intersecting over a line. 
As common in applications of the gauge/gravity duality, 
in the large N, large 't Hooft coupling limit the D3 branes 
are replaced by their classical geometry while the D5 is 
considered as a probe, i.e., its effect on the geometry is 
neglected. The resulting geometry is that of a throat, of 
coordinate length R, pulled by the D3s out of the flat 
10-dimensional (lOd) spacetime: 
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t and x = (x 1 , x 2 , x 3 ) are coordinates on the D3 worldvol- 
ume, A and </> are polar coordinates in the (x 8 , x 9 ) plane, 
p 2 = (x 4 ) 2 + . . . + (a; 7 ) 2 , and dfl 2 is metric on a unit 3- 
sphere. In what follows, we use dimensionless units in 
which R = 1. The D3s are located at A = p = 0, where 
the metric ((2]) has a degenerate horizon. This metric is 
suitable for calculations at zero temperature, the only 
case for which we present detailed calculations here. The 
results, however, do have a bearing on behavior at finite 
T, as we discuss towards the end. 

Matter in the fundamental representation of SU(N) 
(in our case, the electrons) is described by strings stretch- 
ing between the D3s and a probe brane Q. We consider 
the case when the probe D5 wraps x , p and the 3-spherc. 
Since x = x 1 is the only spatial direction shared by the 
D5 and D3s, this brane intersection describes a theory 
of electrons that live on a (l+l)-dimensional defect but 
interact via a (3+l)-dimensional non-abelian gauge field. 
Overall the setup is similar to the system of D3 and D7 
branes intersecting over a plane, in which the electrons 
are confined to move in two spatial dimensions [j| Hoj |. 
Note that in our case there are two directions, x 8 and x 9 , 
orthogonal to all branes. In complex notation, the dis- 
placement of the D5 relative to the D3s in the (x 8 ,x 9 ) 
plane is Ae 1 ^ and forms an order parameter suitable for 
description of superconductivity. The minimal distance 
between the D5 and D3s is the quasiparticle gap (in string 
units). 



For a general D5 embedding, 

A = A(t,x,p,aii) 
(j) = 4>(t,x,p,ai) 



(4) 
(5) 



where on are angles on the 3-sphere. Vibrations of the 
brane correspond to collective modes of the electron fluid. 
In the 't Hooft limit (N — > oo at fixed A = g s N, where 
g s is the closed string coupling) , quantum fluctuations of 
the brane are suppressed by the large value of the brane 
tension, so to the leading order the brane can be consid- 
ered as a classical object. One can then explore various 
embedding ansatzes, which will typically have less coor- 
dinate dependence than the most general form (j4}-([5]). 



We assume throughout that x 
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0. In addition, 



all embeddings we consider here are independent of on. 
With this restriction, the Dirac-Born-Infeld (DBI) action 
of the D5 brane is 

Sdbi - -2ir 2 T 5 J dtdxdpp 3 f/ 4 [-det(G ab + F ab )}^ 2 , 

(6) 

where T5 it the brane tension, G ab for a, b = i, x, p 
are the components of the induced metric, and F ab = 
d a A b — dbA a is a U(l) gauge field (distinct from the usual 
electromagnetic field) on the D5 worldvolume. The clas- 
sical dynamics of the brane is described by the Eulcr- 
Lagrange (EL) equations following from the total of ([6]) 
and a Wess-Zumino term that describes the coup ling of 
the D5 to the background Ramond-Ramond field [111 ]. 

A spatially uniform winding of the order parame- 
ter can be described by (static) embeddings for which 
4>(x) — gx, where q is a constant, while A depends on p 
only. The constant d x <j> sources worldvolume electric field 
F tp through the Wess-Zumino coupling. This leads to a 
system of coupled equations for A and F tp . We do not 
consider this case further here and move on to description 
of a transition at fixed current. 

We begin with static x-independent embeddings A = 
A(p), (j> — 0, At — A t (p) with all other components of 
A a equal to zero. As explained in (in the context of 
a different D-brane system), the value p = A t (oo) is the 
chemical potential for the density ip'ip. Turning to |T]), 
we see that the k > and k < electrons are oppositely 
charged with respect to p. Thus, in the superconduc- 
tor, p is conjugate to the electric current. The lines of 



the radial field F tp have nowhere on the brane to end, 
so the D5 must extend behind the horizon 12j. Hence 
the boundary condition A(p = 0) = 0: at any nonzero 
current the superconductor is gapless. 

The action ([6]) for this type of embedding is Sdbi = 
— J dtJ-, where J- is the free energy: 



T = 2n 2 T 5 J dxdpp 3 ^/f(l + A 2 ( 



(7) 



We use notation A , 



d p A. The Wess-Zumino term 



vanishes. The equation of motion for A t shows that the 
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FIG. 1: Equation-of-state curves for several values of the cur- 
rent J. b is the asymptotic value of the D5 embedding, and 
A iP (0) is its slope at the horizon. The current increases in 
increments of 0.05 from bottom to top. 

current 

J= 2^T 5 SFZ (8) 

is p-independent. Following the usual procedure of Leg- 
endre transforming to go from fixed p, to fixed J, we 
obtain the Legendre transformed free energy 

f = 2n 2 T 5 J dxdp(p e f + J 2 ) 1 /2(l + A^) 1 ^ . (9 ) 

For each value of J, solutions to the corresponding EL 
equation form a one-parametric family with A iP (0) as a 
parameter. The general behavior of a solution at large p 
is a constant: A(oo) = b. 

We refer to the dependence of b on A iP (0) as an 
equation-of-state (EOS) curve. These curves provide a 
convenient way to study the phase transition. Some rep- 
resentative ones are shown in Fig. [T] A nonzero asymp- 
totic value 6^0 corresponds to an explicit breaking of 
the symmetry that shifts by a constant. Since we are 
interested in solutions that break this symmetry sponta- 
neously, rather than explicitly, we look for nontrivial so- 
lutions with 6 = 0. In Fig. [TJ such a solution corresponds 
to a second zero (if any) of the EOS curve; the first zero, 
at A )P (0) = 0, is the trivial solution. We see that sec- 
ond zeroes exist at smaller currents but as J increases 
they decrease in magnitude and eventually merge with 
the trivial solution and disappear. This is the behavior 
characteristic of a continuous phase transition. Numeri- 
cally, the critical value of J is J c = 0.3197. 

To study fluctuations near the critical point, we con- 
sider the linearized EL equations with time and space 
dependence; these can be obtained by Legendre trans- 
forming the quadratic part of the DBI action ^ . Upon 



substitution A(t,x,p) — e 4W *+ lfca: A(p), the linearized 
equation for A reads 

-^a p (VCA, p )+^/oA-^A + ^ = 0, (10) 

where C = p 6 f + J 2 and f = 1 + l/p 4 . At p oo, 
this reduces to the spherical wave equation in (3+1) di- 
mensions. At p — > 0, and u^0, the leading asymptotics 
are A(p) ~ pe ±lu ^ p . We choose the positive sign, corre- 
sponding to waves falling into the horizon. This choice 
will be the source of irreversibility (dissipation) . 

We consider real k and complex lu. Unstable modes 
of the trivial A = embedding are eigenmodes of (fT0|) 
with Im lu > 0. They decay exponentially at both p — > 
and p — > oo. For these boundary conditions, lu 2 is purely 
real; hence lu is purely imaginary. For the moment, let 
us concentrate on the spatially uniform case k 2 = 0. By 
numerically solving (|10|) with Dirichlet boundary condi- 
tions, we find that, when J < J c but close to it, there is 
only one unstable mode, and its profile closely matches 
that of the nontrivial static solution to the nonlinear 
problem (|9|). This is evidence that the instability of the 
trivial embedding develops into one of the nontrivial em- 
beddings we considered earlier. At J = J c , the frequency 
crosses zero in a smooth, analytic manner, approximately 
as Im io — — 0.9 (J — J c ), and the instability disappears. 
Analytically continuing Im lu to J > J c , where the triv- 
ial embedding is stable, we interpret |Im lu\ as a rate at 
which a small displacement of the D5 brane relaxes back 
to A = 0, i.e., a superconducting fluctuation decays into 
normal electrons. 

Scaling of the relaxation rate with the control parame- 
ter near a critical point is described by a dynamical crit- 
ical exponent, usually called y. The linear dependence 
of the rate on J — J c found above corresponds to y = 1 . 
This, of course, is not the true infrared scaling. The lat- 
ter is determined by interaction among the modes of (|10j) : 
however weak initially, this interaction may renormalize 
to large strengths at large distances. The classical ap- 
proximation to the D5's dynamics (corresponding to the 
leading large N approximation on the gauge theory side) 
does not see this renormalization but it can be used to 
determine the strength of the interaction — and, impor- 
tantly, its sign — at the cutoff scale. To do that, we first 
consider a static but x-dependent configuration 

Ae i<t '(x,p; J) = J dke lkx ^(k)A^{p;k, J) , (11) 

where ^ is a complex amplitude, and A( ) is the unstable 
mode of (fT0)l at the corresponding values of the parame- 
ters (J < J c ). Substituting this into the DBI action (j6]) 
and expanding to the fourth order in '3/ (but only to the 
second for terms containing derivatives with respect to 
x) , we obtain a Ginzburg-Landau-type free energy. With 
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a suitable normalization of "J, it reads 



dkn Q (k, J)|*(fc)| 2 + c dx\*(x)\ 4 



/(12) 

where S7o is Im lo of the unstable mode, and *£(x) is the 
Fourier transform of *&(&). Slo ~ a (^c — </) — a'fc 2 for J 
near J c and small fc 2 , while c is approximately constant. 
The constants a, a', and c are all positive. 

The rate of instability growth (at J < J c ) or relaxation 
(upon continuation to J > J c ) can be encoded into an 
effective action by allowing a weak time dependence in 
'F, going over to the Euclidean time via t = — ir, and 
adding the dissipative term 



issip 



7tT 5 / dxdn\n\\t>(x,n)\' 2 



(13) 



where f2 is the Euclidean frequency. The full Euclidean 
effective action is Se = J drToh + Sdissip- This theory 
is known as the dissipative XY model. It has been con- 
sidered in applications to Id superconductors [13| but, as 
far as we know, not specifically in reference to transition 
at a critical current. The gravity dual in effect provides 
a microscopic derivation of this model for that particular 
case. 

As seen from Q12p. (IT31 . the large T5 (in our dimen- 
sionless units, T5 ~ yXN) suppresses fluctuations of vF 
Thus, the ultraviolet behavior of the system is controlled 
by the Gaussian fixed point c = 0. It is known that 
this fixed point is unstable under renormalization, and 
the theory flows to a nontrivial infrared fixed point (for 
a summary of its properties, see Ref. [Hj]). However, 
since the true critical behavior will require a large length 
scale to set in, even a small deviation of J from J c may 
lock the system in the domain controlled by the Gaussian 
fixed point. That is the case when the renormalization, 
say, for J < J c is cut off at the finite (correlation length) 
scale £ ~ £o(l — J/Jc)^ 1 ^ 2 , where £0 is some microscopic 
length, while the quartic coupling remains small. 

Experimentally, superconducting nanowires have been 
observed to switch to the normal state at currents below 
the estimated critical current 0-0] • These experiments 
have found large fluctuations in the switching current, 
which persist down to low temperatures and which have 
been interpreted as a consequence of phase slips. 

This interpretation is consistent with our conclusion that 
the transition at J = J c is second-order, as it implies 
that the free energy barrier suppressing phase slips can 
be almost completely removed by bringing J close to J c . 
This in turn implies that the transition at J — J c must 
be continuous or at most very weakly first-order. 



The following estimate suggests that the near-critical 
behavior in these experiments is controlled by the unsta- 
ble Gaussian fixed point, rather than the stable XY one, 
at least at temperatures where thermally activated phase 
slips are thought to be the main effect. The exponential 
factor in the rate of thermal activation is exp(— SF/T), 
where 5 J 7 is the free energy barrier and T is the tem- 
perature. Near J = J c , S T scales as the product of 
the free energy density and the correlation length, i.e., 
as ( J c — J) 2 ~ a ~ v , in terms of the conventionally defined 
critical exponents. For the Gaussian point, a = and 
v = |, so 5 J 7 ~ (J c — J) 3 / 2 . Curiously, this is the same 
scaling as obtained for a Josephson junction 14| . It has 
been found to provide a good fit to the data in Ref. @ 
and for the amorphous samples in Ref. Q (for the crys- 
talline samples [7|, the 5/4 power law has been found to 
be a better fit). In contrast, for a fixed point obeying 
hyperscaling, 2 — a — v = and SJ 7 scales to a constant. 
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